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ABSTRACT 

For  the  numerical  computation  of  efficiencies  for  optical  gratings,  there  exists 
a huge  variety  of  algorithms.  Dealing  with  a boundary  value  problem  for  an 
elliptic  partial  differential  equation,  the  application  of  finite  element  methods 
(FEM)  is  natural  too.  However,  the  oscillatory  nature  of  the  electromagnetic 
fields  requires  some  modifications.  The  resulting  FEM  program  can  be  used  as  a 
part  of  an  algorithm  to  design  optimal  gratings. 

THE  FINITE  ELEMENT  METHOD  (FEM) 

The  variational  form  of  the  boundary  problems  is  well  known  and  its  coercive- 
ness is  well  established  (cf.  e.g.  [6,1,3]).  For  example  the  variational  equation 
for  u(x,y ) = v(x,y ) ■ exp(-iax)  with  v the  unknown  third  component  of  the 
amplitude  of  the  scattered  magnetic  field  in  the  case  of  TM  polarization  is 

J ^{V  + *(«>  °)}«  ' (v  + i(«>  °)}<F  - ulP+  ^+0C' 

+yhl-(T°u)* = -?Fp  jjM***n-  *> e W « 

The  domain  Q is  the  rectangular  cross  section  of  the  grating  profile  taken  over 
one  period  and  T*  stand  for  the  upper  resp.  lower  boundary  sides.  The  symbol 
k stands  for  the  piecewise  constant  refractive  index  taking  the  constant  values 
k+  resp.  k~  above  resp.  below  the  grating.  The  number  a is  the  product  of  k+ 
times  the  sine  of  the  incidence  angle,  V is  the  gradient,  and  a hypersingular 
boundary  integral  operator.  The  test  function  c p runs  through  all  periodic  (w.r.t. 
to  the  first  variable  x)  functions  in  the  Sobolev  space  JHrl(f2).  A variational 
equation  similar  to  (1)  holds  for  the  TE  polarization.  For  the  case  of  conical 
diffraction  (incident  plane  wave  with  direction  not  in  the  plane  perpendicular  to 
the  grooves  of  the  grating) , a coupled  system  of  two  such  variational  equations  is 
to  be  solved. 

To  get  the  FEM  solution,  the  domain  is  split  into  triangles/rectangles.  An 
approximation  of  the  magnetic  field  is  sought  in  form  of  a continuous  piecewise 
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linear/bilinear  function,  and  the  test  functions  in  (1)  are  replaced  by  the  piece- 
wise  linears/bilinears.  Substituting  the  numerical  solutions  for  u over  T*  into 
well-known  integral  representations,  we  get  the  reflected  and  transmitted  energy 
and  the  efficiencies  of  the  reflected  and  transmitted  modes.  In  contrast  to  pop- 
ular alternative  methods,  the  FEM  needs  no  Rayleigh  expansion,  no  slicing  of 
blaze  gratings,  and  no  solution  of  big  dense  systems  of  integral  equations.  The 
convergence  of  the  method  is  well  established. 


Example  : cover  material  : 

coating  material  : 
substrate  material  : 
polarization  : 
period  : 

incidence  angle  : 
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The  figures  show  a coarse  FEM  grid  of  a coated  lamellar  grating  profile  with 
overetching,  the  isolines  of  the  real  part  of  the  solution  u,  and  the  efficiencies  of 
the  transmitted  modes  of  order  -2,  -1,  0,  and  1 depending  on  the  incidence  angle 
in  the  interval  [40°,  50°]. 


GENERALIZED  FEM 

Unfortunately,  due  to  a mismatch  of  the  frequency  of  the  solution  and  the  discrete 
frequency  of  the  approximate  solution,  the  FEM  deteriorates  with  large  wave 
numbers.  Even  if  the  function  can  be  approximated  in  the  space  of  piecewise 
linears  with  high  accuracy,  the  error  of  the  FEM  solution  of  Eq.  (1)  may  be 
large.  To  overcome  this  so-called  “pollution  effect”,  we  have  first  implemented 
a generalized  FEM  for  the  case  of  lamellar  (binary)  gratings  which  is  a finite 
difference  scheme  over  uniform  rectangular  gratings  (cf.  [2,4]).  Further,  we  have 
tested  the  partition  of  unity  method  together  with  mortar  techniques  (cf.  [5]), 
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and  we  have  implemented  a new  generalized  FEM  based  on  trial  functions  the 
restrictions  of  which  to  the  triangles  of  the  triangulations  are  solutions  of  the 
partial  differential  equation.  These  solutions  are  generated  by  an  FEM  over  a 
fine  uniform  grid  in  the  subtriangle.  For  a fixed  accuracy,  these  methods  reduce 
the  storage  requirements  and,  in  some  cases,  the  computing  times  essentially.  In 
the  table  we  present  results  for  a coated  echelle  grating  (wave  length  160  nm,  166 
grooves  per  mm,  width  of  MgF2  coating  25  nm,  blaze  angle  80°,  apex  angle  90°, 
cover  material  Air,  grating  material  Al,  TE  polarization,  incidence  angle  80°). 


degrees  of  freedom 

memory  for  solver 

efficiency  of  order  74  (refl.) 

105  785 

0.35  GB 

37.931045 

263624 

0.70  GB 

67.384460 

559800 

1.98  GB 

68.390312 

OPTIMZATION  OF  GRATINGS  (SYNTHESIS  PROBLEM) 

The  mentioned  methods  for  the  numerical  solution  of  the  direct  diffraction  prob- 
lem can  be  used  as  a part  of  an  algorithm  to  design  optimal  gratings.  We  have 
implemented  a code  (cf.  [4])  to  minimize  several  objective  functions  including  ef- 
ficiencies and  phase  shifts.  On  a set  of  coated  lamellar  grids  containing  a certain 
number  of  rectangular  pieces  with  prescribed  material  properties,  we  determine 
an  optimal  grating  by  a gradient  descent  method.  The  latter  is  based  on  the 
efficient  computation  of  the  gradients  by  generalized  FEM. 

ACKNOWLEDGEMENT 

The  authors  gratefully  acknowledge  the  support  of  the  German  Ministry  of  Ed- 
ucation, Research  and  Technology  under  Grant  No.  03-ELM3B5. 

REFERENCES 

[1]  G.  Bao,  D.C.  Dobson,  and  J.A.  Cox,  “Mathematical  studies  in  rigorous  grat- 
ing theory”,  J.Opt.Soc.Amer.  A 12,  1029-1042  (1995). 

[2]  I.  Babuska,  F.  Ihlenburg,  E.  Paik,  and  S.  Sauter,  “A  generalized  finite  element 
method  for  solving  the  Helmholtz  equation  in  two  dimensions  with  minimal 
pollution”,  Comp. Methods  Appl.Mech.Eng.  128  325-359  (1995). 

[3]  J.  Elschner,  R.  Hinder,  F.  Penzel,  and  G.  Schmidt,  “Existence,  uniqueness 
and  regularity  for  solutions  of  the  conical  diffraction  problem” , Math. Models 
& Methods  Appl.Sci.  10,  317-341  (2000). 

[4]  J.  Elschner  and  G.  Schmidt,  “Numerical  solution  of  optimal  design  problems 
for  binary  gratings",  J.Comput.Phys.  146,  603-626  (1998). 

[5]  J.M.  Melenk  and  I.  Babuska,  “The  partition  of  unity  method:  Basic  theory 
and  applications”,  Comput. Methods  Appl.Mech.Eng.  139,  289-314  (1996). 

[6]  H.P.  Urbach,  “Convergence  of  the  Galerkin  method  for  two-dimensional  electro- 
magnetic problems”,  SIAM  J.Numer.Anal.  28,  697-710  (1991). 


Kiev,  Ukraine,  IX-th  International  Conference  on  Mathematical  Methods  in  Electromagnetic  Theory 


